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We discuss random matrix models in terms of elementary operations on Blue's 
functions (functional inverse of Green's functions). We show that such operations 
embody the essence of a number of physical phenomena whether at /or away from 
the critical points. We illustrate these assertions by borrowing on a number of 
recent results in effective QCD in vacuum and matter. We provide simple physical 
arguments in favor of the universality of the continuum QCD spectral oscillations, 
whether at zero virtuality, in the bulk of the spectrum or at the chiral critical 
points. We also discuss effective quantum systems of disorder with strong or weak 
dissipation (Hatano-Nelson localization). 
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1. INTRODUCTION 



Complex and disordered systems play an important role in a number 
of areas in physics ranging from mesoscopic conductance fluctuations to 
simplicial quantum gravity. The dual interplay between disorder and lo- 
calization |U], and disorder and chaos Q, has far reaching consequences on 
many aspects of modern physics whether at the macroscopic or microscopic 
level. 

An important tool for investigating Hamiltonians of disordered systems 
is random matrix theory ||, ||, [| . A number of techniques have been de- 
veloped such as the method of orthogonal polynomials, the supersymmetric 
and replica methods to cite afew. A powerful alternative to some of these 
methods has been the work of Voiculescu in the context of operator alge- 
bras and the ensuing concept of free random variables. Free disorder allows 
for an exact calculation of various moments of a Hamiltonian composed of 
random plus deterministic parts, hence it's Green's function, while main- 
taining its parts logically separated. In physical systems, freeness amounts 
to summing over all 'single-site' rescatterings ||, a procedure that is close 
to mean-field treatments in many-body physics. In systems amenable to a 
1 /N expansion, this procedure is equivalent to the one based on the concept 
of independent disorder in the large N limit. In most cases, it is simpler to 
implement and holds for arbitrary forms of disorder. 

In these lectures we will review and generalize the basic operations in- 
troduced by Voiculescu for free random variables, and show how they could 
be used in a variety of physical applications. Instead of using the standard 
language of Green's functions, we choose to work with their functional in- 
verses introduced and named Blue's functions by Zee ||. We show that the 
use of Blue's functions sheds new light on a number of issues in conventional 
as well as non-conventional random matrix theory, while providing an effec- 
tive calculational shortcut. In many ways, it is the user- friendly alternative 
to Voiculescu's involved concepts of non-commutative independence (jjj and 
non-crossing cumulants to cite a few The concept of addition law be- 
tween Blue's functions in random matrix theory is perhaps the closest in 
spirit to the concept of a mean-field analysis in many-body physics. In- 
deed, for quantum problems with diagonal disorder it is equivalent to the 
Coherent Phase Approximation (see below and |lT|| ). 

In the first part of these lectures we will introduce the elementary cal- 
culus with Blue's functions. We will define the concepts of addition law and 
multiplication law for Blue's functions for hermitian random ensembles with 
arbitrary measures, and the concept of generalized addition law for gener- 
alized Blue's functions for nonhermitian random ensembles with arbitrary 
measure. We also define the concepts of differentiation and integration of 
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Blue's functions and their respective relationship with the end-points and 
boundaries of eigenvalue spectra. Finally we exploit some analytical prop- 
erties of Blue's functions by constructing conformal mappings between the 
domains of analyticity for hermitian and non-hermitian ensembles. In the 
second part of these lectures, we will present a number of applications to 
effective Quantum Chromodynamics (QCD), and models of quantum me- 
chanics with strong and weak dissipation. 

In effective QCD we emphasize the role played by the quark modes near 
zero virtuality, in the characterization of the issue of spontaneous symmetry 
breaking in the vacuum and near the chiral critical points. We present 
simple arguments for the importance of chiral symmetry in the conditioning 
of the spectral fluctuations near zero virtuality as well as within the bulk 
of the low-lying eigenvalue spectrum, both in vacuum and matter. We 
introduce simple examples of chiral random matrix models, and provide 
plausible arguments for the decoupling of the soft and hard modes in infrared 
sensitive observables. We give a schematic description of the U(l) problem 
and its resolution, and illustrate the ingredients at work in the screening 
of the topological charge. We also show how the concepts of temperature 
and chemical potential may be schematically incorporated in the present 
approach, and how the structural changes in the eigenvalue spectra are 
encoded in the Blue's functions. We briefly comment on the role of masses 
and vacuum angle 9 in this effective approach to QCD. In the models of 
quantum mechanics, we demonstrate the versatility of our approach to the 
problem of open chaotic scattering as well as tunneling in the presence of 
dissipation. We go on showing that the addition law for Blue's functions for 
a class of quantum systems with dissipation is the analogue of the Coherent 
Phase Approximation (CPA), a method commonly used for the Anderson 
model Q and its variants O, We speculate on further applications of 
this approach in our conclusions. 



2. OPERATIONS on BLUE'S FUNCTIONS 



The fundamental problem in random matrix theories [|, g, 13, |14| is to 
find the distribution of eigenvalues in the large iV (size of the matrix A4) 
limit. One encounters two generic situations — in the case of hermitian 
matrices the eigenvalues lie on one or more intervals on the real axis, while 
for general non-hermitian ones, the eigenvalues occupy two-dimensional do- 
mains (or even fractals) in the complex plane. In both cases the distribution 
of eigenvalues can be reconstructed from the knowledge of the Green's func- 
tion 
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where averaging is done over the ensemble of N x N random matrices gen- 
erated with probability 

P( M ) = i e -"*v(M) _ (2) 

For hermitian matrices the discontinuities of the Green's function coincide 
with the support of eigenvalues. In the simplest example of the Gaussian 
potential V{M) = \M. 2 for a hermitian matrix M, a standard calculation 
gives @ 

G{z) = \{z-yf#=l). (3) 

The reconstruction of the spectral function (eigenvalue distribution) is based 
on the well known relation 

— — - = PV- T iir5(x) (4) 
x ± le x 

hence 

v(\) = -- \m\ lm.G{z)\ z= x +it . (5) 

In the Gaussian case (|3|), the discontinuities in the Green's function come 
from the cut of the square root, leading to Wigner's semicircle law 2ttv{X) = 
y/4 — A 2 for the eigenvalue distribution. In practice, the averaging procedure 
< ... > for more complicated probability distributions is difficult. It could 
be simplified if the randomness is assumed to be free (see below). 

In the general case of non- hermitian matrices, the two-dimensional do- 
mains where the eigenvalues are distributed are exactly those where the 
Green's functions cease to be holomorphic, that is 

v(z,z) = -dgG(z,z)£0. (6) 

TT 



2.1. Addition 



Let us consider the problem of calculating Green's function for the sum 
of two independent ensembles M\ and M.2, i.e. 

G(z) = 1 J [dM 1 ][dM 2 }Pi(M 1 )P 2 (M 2 ) tv z _ M \_ M2 
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The concept of addition law for hermitian ensembles was introduced by 
Voiculescu [0. In brief, he proposed an additive transformation (R trans- 
formation), which linearizes the convolution of non-commutative matrices 
([?]), much like the logarithm of the Fourier transform for the convolution 
of arbitrary functions. This method is an important shortcut to obtain the 
equations for the Green's functions for a sum of matrices, starting from the 
knowledge of the Green's functions of individual ensembles of matrices. This 
formalism was reinterpreted diagrammatically by Zee |]], who introduced 
the concept of Blue's function. 

Let us consider the problem of finding the Green's function of a sum 
of two independent (free 0) random matrices Mi and M2, provided we 
know the Green's functions of each of them separately. The key observation 
is that the self-energy S, defined as £(2) = z — G~ 1 (z), can be always 
expressed as a function of G itself and not z as usually done in textbooks. 
Then the following addition law holds 

Ei+ 2 (G) = Si(G) + Si(G). (8) 

This formula does not hold for G — > z. Voiculescu's R transform is nothing 
but the self-energy expressed as a function of G {i.e. R(G) = E(G)). For 
the Gaussian randomness R(u) = u. For an arbitrary complex number u, 
the addition (|8|) reads R% + 2(,u) = Ri(u) + R2(u). The set of all R transforms 
generate an Abelian group. The Blue's function, introduced by Zee ||, is 
simply 

B(G) = S(G) + G~ x . (9) 

Therefore, using the identity G(z) = [z — X)" 1 , we see that the Blue's 
function is the functional inverse of the Green's function 



B[G(z)) = z 
and the addition law for Blue's functions reads 



(10) 



B 1+2 (z) = B x {z) + B 2 {z) 



(11) 



Using the definition of the Blue's function we could rewrite the last equation 
in "operational" form as 



Bi(G) + B 2 (G) - G- 1 



(12) 



with G1+2 = G. The algorithm of addition is now surprisingly simple ||: 
Knowing Gi and G2, we find Bi and B2 from (10), and read from (12) the 
final equation for the resolvent of the sum. Note that the method treats on 
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equal footing the Gaussian and non-Gaussian ensembles, provided that the 
measures Pi and P2 are independent (free). 

The generalization of this algorithm to the case of arbitrary complex 
random matrices is subtle. In this case, the eigenvalues are complex and 
are not distributed along the one dimensional cuts on the real axes, but 
rather form on the z plane two-dimensional surfaces or fractals fll7jj . Again 
we would like to have a method for finding the Green's function (|l|). 

The generalization |l8|] amounts to consider, in addition to the original 
matrix M., a hermitian conjugate copy and an infinitesimal 'coupling' A 
between both copies. This leads us to investigate the matrix-valued Green's 
function 



Gqq Gqq \ / f Z - M A 



Gqq Gqq ) 2 Nx2N \V A Z / 2.\ -<2.\ I 



(13) 



The subscript 2N x 2N indicates that we are dealing with a 2N by 2N 
matrix — M. and are N hy N matrices, and z and A are implicitly 
multiplied by an N x N identity matrix. A variant of this construction 
leading to similar results was recently proposed in [19]. Now, the ordinary 



Green's function (|T]), the chief aim of the calculation, is just the trace of the 
upper left-hand corner of G- 

G(Z) = lim ^tl N Gqq (14) 

JV— too iV 

where the limit N — > 00 is understood before A — > 0. However in ordinary 
calculations one never has to perform any explicit limits. It is convenient 
to take the traces of the four N by N blocks of G to obtain a 2 x 2 matrix: 



Gqq Gqq \ ( tr^f Gqq trjv Gqq 



Gqq G m J ^ \ tl N Gqq tVN Gqq J 2x2 



(15) 



The generalized Blue's function ffjq , 20 1 is now a matrix valued function 
of a 2 x 2 matrix variable defined by 

B(G)=Z=( Z X J) (16) 

where A will be eventually set to zero. This is equivalent to the definition 
in terms of the self-energy matrix 

B{G) = ^ + Q- 1 (17) 
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where £ is a 2 x 2 self energy matrix expressed as a function of a matrix 
valued Green's function. The generalized addition law reads 

Z = B^G) + B 2 {G) - Q- 1 . (18) 

This matrix equation has always two kinds of solutions. The first one (triv- 
ial) corresponds to the case when the off-diagonal matrix elements of Q are 
equal to zero. In this case, (|T^ ) splits into two copies of the addition law 
(|IT|), for the holomorphic function B(z) and its anti-holomorphic copy B(z), 
respectively. The second (non-trivial) solution corresponds to the case when 
holomorphic separability of ( |18| ) is no longer possible, and all matrix ele- 
ments of Q are nonholomorphic functions of z and z. The average density of 
eigenvalues follows then from the electrostatic analogy to two-dimensional 
Gauss law [16| 



u(z, z) = -d 2 G qq (z, z) . (19) 

The condition that the holomorphic and nonholomorphic solutions match 
(i.e. G qq -(z, z) = 0) defines the borderline of the eigenvalue distribution of 
the z plane. The advantage of the addition law, for hermitian and non- 
hermitian random matrix models, stems from the fact that it treats Gaus- 



sian and non-Gaussian randomness on the same footing |21]. 



2.2. Multiplication 

We now demonstrate that Blue's functions provide also an important 
shortcut to obtain the equation for the Green's (Blue's) function for a prod- 
uct of matrices, starting from the knowledge of Green's (Blue's) functions 
of individual ensembles of matrices, i.e. to find 

G(z) = i-/tr ~ ~ — \ (20) 

w N\ z-Mi -Mil 

provided that .Mi and A4 2 are free and (trA^i), (trA^) ^ 0. Let us intro- 
duce the notation for continued fraction 



w w 



(21) 



Then the multiplication law for Blue's functions reads 
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For the diagrammatic proof of this relation we refer to [22]. Contrary to 
standard Green's functions, the Blue's functions provides the factorization 
mechanism for the averaging procedure in (20). Using the definition B{G) = 
Y,(G) + G~ 1 we get B{w /£(•)) = £(»)(l + l/u>) therefore the multiplication 
law reads 

Note that the S-transformation of Voiculescu is related to £(•) via 22] 

S( W ) = -1 (24) 

hence (|2^) reads Si^iw) = Si(w) ■ S2(w). 

To check that the relation (]24}) with Voiculescu's S-transform is true, 
we have to check that the function x( w ) defined by S(w) = x(w;)(l + 1/w) 
satisfies J?]] 

x(zG(z)-l) = l/z. (25) 
Comparing the definition of x an( i (13) we see th & t indeed 



1 1 1 1 

x[w) ~ (i + i/^s^) " %M " s[gRJ = " (L,,,) 



z 



where in the last two equalities we changed the variables w = zG(z) — 1 and 
used the definition of Blue's function, respectively. 



2.3. Differentiation 

In the previous two sections we showed how to use Blue's functions to 
obtain the Green's function of a sum or product of random matrices from 
the properties of the individual ensembles. Now we will concentrate on the 
information about the structure of the eigenvalues of some fixed random 
matrix ensemble, which we may extract from the Blue's function. 

As shown by Zee ||, the differentiation of the Blue's function leads to 
the determination of the endpoints of the eigenvalue distribution for her- 
mitian random matrices. One could easily visualize this point recalling the 
form of the Green's function for the Gaussian distribution (|3|). When ap- 
proaching the endpoints z = ±2, the imaginary part of the Green's function 
vanishes like >/(2 — A)(2 + A). Therefore, the endpoints fulfill the equation 
G' (z)\ z= ±2 = oo, which may be used as the defining equation for unknown 
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endpoints in the case of more complicated ensembles, since the endpoints 
are the branching points of the resolvent. Since Blue's function is the func- 
tional inverse of the Green's function, the location of the endpoints could 
be inferred from its extrema 

dB(G) 

-&■ = <>. (27) 



2.4- Integration 

In this section we will consider the phase structure of a class of systems 
with a partition function depending on a parameter z: 

Z N (z) = (det(>- M)). (28) 

In general A4 may be a sum of a random matrix and some deterministic 
one (not necessarily hermitian). The variable z may also be promoted to 
a matrix. A typical example being the chiral random matrix model for 
effective QCD as discussed below. Generically, the analytical structure of 
(f28|) and its zeros are of interest for phase studies in generalized statistical 



mechanics [23]. In matrix models, this information is all contained in the 
Blue's functions. 

For finite ./V the partition function is a polynomial in z and has N (Yang- 
Lee) zeroes 

Z N (z) = (z - z\)(z - z 2 ) ■ ■ ■ (z - z N ) . (29) 

Taking the logarithm and approximating the density of zeroes by a contin- 
uous distribution we get 

log <det(z -M))=Y J log(* -Zi) = J p(z') log(z - z')dz' . (30) 

i 

We see that the density of zeroes can be reconstructed from the discontinu- 
ities of the unquenched Green's function 

d z log (det(z - M)) = (tr — V — ■ det(z - M)\ (31) 

through Gauss law [24|. For z close to infinity the unquenched and quenched 
(||) Green's functions coincide^. Moreover the unquenched resolvent is non- 
singular configuration by configuration, hence holomorphic. Therefore it is 
exactly the functional inverse of the ordinary holomorphic Blue's function. 

1 In leading order in 1/N. 
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In general we may have several branches for the Green's functions, which 
are determined by different saddle points contributing to (]28|) . The discon- 
tinuity (cusp) and hence the location of the zeroes is determined by the 
possible contribution of two of them. Prom a 1/N expansion 

\ogZ N = NE + E 1 + ^E 2 + ... (32) 

we see that two saddle points may contribute if 

ReE s Q pJ = ReE s Q pJI (33) 
and Eq is determined by ( |3l| ) 

E = J dz'G(z') + const (34) 



or equivalently 



E = zG- JdG B(G) + const (35) 



after integrating by part. Note that we have used the fact that z(G) = B{G) 
is just the Blue's function of G. The constant in Eq is fixed by the asymp- 
totic behavior of (p3), that is Zjy ~ z . We should add that performing 
the integration using Blue's function is in most cases much simpler than the 
corresponding direct integration of G(z). 

If we denote the complex valued set of all the cusps by 

F(z, z) = ReE s Q pJ - ReE s pJI = (36) 

then, to leading order, the distribution of singularities along the "cusps" 
(H is 

6 (z) = ±-\d z F\ 2 6(F(z,z)) (37) 

which is normalized to 1 in the z-plane. Redefining the density of singu- 
larities by unit length along the curve F(z, z) = 0, we may rewrite ( [37]) 
as 

6 (z)\ F=0 = ±\d z F\ = ±-\GU-GM\. (38) 

For q ^ 0, the integrand in ([3(]) is singular at z = z' which results in a cusp. 
For q = 0, that is d z F = 0, (|30|) is differentiable. For physical Eq (real 
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and monotonically increasing), the points g = are multi-critical points. 
At these points all n-point (n > 2) functions diverge. This observation also 
holds for Ising models with complex external parameters [E9]. Assuming 



macroscopic universality [26, 27, 28, 29, [30| for all n-point functions(n > 2), 



we conclude that d z F = means a branching point for the resolvents, hence 
d z G = co or B'(G) = 0, as shown before. 

The next to leading contribution E\ in ( |32"|) does not influence the pre- 
ceding discussion except at exceptional pointsg. In fact the expression for 

Zjq is 

Z N (z) = e NEo -({detC(G)}^+o(^ . (39) 

The explicit form of the argument of the determinant could be addressed 
using several methods. We do not discuss this point here. We would like to 
stress, that usually C is a simple function of the resolvent G on a specific 
branch. In such a way the Blue's function conditions the leading behavior 
of the partition functions, either explicitly via integration leading to Eq, 
or implicitly via the functional inverse and the macroscopic universality of 
higher correlators. 

Note that due to the power —1/2 in (|39|), the fluctuations have "bosonic" 
character, and are dwarfed by the "quark-loop" contribution in the ratio 



1 : N [31]. In particular, the condition for the divergence of the "bosonic" 
fluctuation at C(G) = 0, signals the breakdown of the 1/N expansion. These 
are the exceptional points, for which new regimes of microscopic universality 
may set in, following pertinent scalings. The standard example for Gaussian 
randomness is the divergence of the wide two-point correlator for z 2 — 4 = 0, 
signaling the universal behavior (Airy oscillations) |32| near the endpoints 
of the spectra. In the case of chiral Gaussian randomness the divergence 
happens for z(z 2 —4) =0. The additional multicritical point z = is due to 
the chiral nature ("Goldstone pole") and signals another class of universal 
behavior (Bessel oscillations) 



2.5. Mapping 

The concept of Blue's functions has a very surprising interpretation in 
terms of conformal mappings between two distinct ensembles of random 
matrices. Indeed, these mappings can take the discontinuities or cuts of the 
eigenvalue distribution of a hermitian random matrix model and transform 
them into the boundary of the eigenvalue distribution of its non-hermitian 

2 For a detailed discussion of this point see [IL8|. 
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analogue |20| . These transformations can be carried fully in the holomorphic 
region, thereby avoiding the subtle issue of nonholomorphy altogether. 

To illustrate these points, consider the case where a Gaussian random 
and hermitian matrix H is added to an arbitrary deterministic matrix M. 
The addition law for Blue's functions states that 

B h+ m (u) = B h (u)+B m (u)--=u + B m (u) (40) 

it 

where we have used explicitly that for a Gaussian ensemble Bjj(u) = u + 
1/u. Now, if we were to note that in the holomorphic domain the Blue's 
transformation for the Gaussian anti-hermitian ensemble is Bih(u) = — u + 
1/u, we fmd0 

B iH+M (u) = B iH (u) + B M (u) - - = -u + B M (u). (41) 

u 

These two equations yield 

B iH+M (u) = B H+M {u) - 2u . (42) 

Substituting u — > Gh+m{z) we can rewrite (fffi) as 

BiH+M[G H +M(z)} = z - 2G h +m(z) . (43) 

Let w be a point in the complex plane for which Gih+m{w) = Gh+m(z). 
Then, using the definition of the Blue's function, we get 

w = z-2G h +m(z). (44) 

The result ( |4"4l ) provides a conformal transformation mapping the holomor- 
phic domain of the ensemble H + M (i.e. the complex plane z minus cuts) 
onto the holomorphic domains of the ensemble iH + M, (i.e. the com- 
plex plane w minus the "islands"), defining in this way the support of the 
eigenvalues. 



3. UNIVERSALITY without BLUE'S FUNCTIONS 

The concept of Blue's functions rely on the use of the 1/N expansion. 
There is a number of circumstances where the standard 1/iV expansion 
breaks down, and where interesting phenomena take place in spectral anal- 
ysis. In this part of our lectures, we will discuss them in the context of 

3 Note that in the holomorphic domain, anti-hermitian Gaussian nullifies the hermitian 
Gaussian in the sense of group property for addition law for the R transformation. 
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Figure 1: Typical Dirac Spectrum for massless chiral quarks in Euclidean 
space (a) . The enlargement around zero virtuality for finite volume is shown 
in (b). 



effective QCD. The Euclidean Dirac spectrum in effective QCD presents 
the global structure described qualitatively in Fig. [l]. For gauge configura- 
tions with non-zero winding number n / 0, an accumulation of zero modes 
take place at zero as expected from the Atiyah-Singer index theorem. The 
large eigenvalues A do not sense sufficiently localized gauge configurations. 
They are described by a free density of states z/(A) ~ ^|A| 3 with V the Eu- 
clidean four-volume. Near zero, the density of states is conditioned by the 
spontaneous breaking of chiral symmetry. This phenomenon causes a huge 
accumulation of eigenvalues, with a level spacing typically of order 1/V. In 
this regime the 1 /N expansion has to be reorganized to accomodate for new 
scaling laws. We now proceed to describe some of these features and other 
related properties, using effective QCD. 



3.1. Effective QCD and Broken Chiral Symmetry 



Of all the effective model approaches to QCD, most noteworthy are those 
that account properly for the spontaneous breaking of chiral symmetry []34| . 
There are a number of effective models in this direction, such as the linear 
and non-linear sigma models and their progenies, the Nambu-Jona-Lasinio 
model and its variants [35|, the instanton model [34, 36], and more recently 
the chiral random matrix models.^] Since these lectures are devoted to ran- 
dom matrix models, we will focus our attention on the latters. 

Why do chiral random matrix models have anything to do with effective 
QCD? The answer to this question is simple: in the long wavelength limit, 
chiral symmetry is spontaneously broken in the vacuum. As a result the 
symmetry of the associated coset space as well as the mode of explicit chiral 

4 We note at this stage that chiral perturbation theory |37]] or chiral reduction formu- 
lae J3fi[] , are not models to QCD but a way to analyze its chiral content using (broken) 
chiral symmetry and data. 
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breaking determines uniquely the character of the QCD effective action in 
the infinite wavelength limit and to leading order in the symmetry breaking. 

For Nf > 1 the coset space is isomorphic to SU(Nf) and the mode of 
explicit breaking is in the fundamental (Nt,Nf) representation (ignoring 
the induced U(l) breaking from the anomaly []). The partition function of 
four dimensional QCD on an Euclidean hypertorus of volume V with finite 
vacuum angle 9 is |[7|, |5(| 

Z(0, M) = ZJ [dU]e~^ ^Me^fU + U^Me-^I )+ 0) (4g) 



where the mass matrix M is real, positive and diagonal, and U is SU(Nf) 
valued. O are the nonzero mode contributions as well as terms of order M 2 
and higher. Zq is an overall M, 8 independent normalization, and V is the 
space-time volume. S = | (qq) | is a measure of the vacuum condensate in 
the chiral limit with a normalization scale matching that of the quark mass 
matrix []. Chiral random matrix models share (|45|) with QCD ^l], |33| . 
(This is also the case of all the other models mentioned above.) 



In the limit VMS 1 (macroscopic regime), the integrand in (45) is 
peaked. As a result, the contribution is dominated by the saddle point 
with a preferred direction on the coset, say U ~ 1. In this regime, chiral 
symmetry is spontaneously broken, with a vacuum condensate of the order 
of S (in absolute value). The saddle point contribution is extensive with a 
vacuum energy density £ = MS + O. The corrections O may be organized 
in chiral power counting My/V ~ 1. The additional tree contribution to 
the saddle point due to the omitted non-zero modes is of order 1/y/V, and 
vanishes in the large volume limit. Similarly, the loop contribution from 
the non-zero modes (massive bosons) as well as the counterterms to order 
1/V vanishes as V — > oo. The former disappears exponentially through 
Boltzman-like factors e~^^^ with /? 4 = V. The divergences in the loops 
are solely extensive (same in small and large volumes) and readily absorbed 
in the renormalization of the cosmological constant. Chiral random matrix 
models support this rationale only to tree level in the macroscopic limit, 
and hence should be regarded as effective models in this regime. 

In the limit VMT, -C 1 (microscopic regime), the condensate is of the 
order of the volume V, and vanishes in the chiral limit (Mermin-Wagner- 
Coleman theorem). The calculation can still be organized in powers of M 
and 1/V such as MV ~ lQ. In this case ( [IB] ) is still the leading contribu- 

5 The inclusion of the U(l) anomaly upsets the chiral power counting, except in large 
N c where the anomaly is treated perturbatively |39| , [to| 

6 For Nf = 1 chiral symmetry is explicitly broken by the U(l) anomaly. Nevertheless, 
power counting in M is still valid and the representation ( jiil ) on the U(l) coset still holds, 
provided that the term Vx*(i\n det U) 2 is added p9| . 

7 This is just a reorganization of the chiral power counting. 
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tion from the zero modes (constant modes on the coset). In this way of 
counting, the non-zero mode contributions are Gaussian and subleading in 
the large volume limit, making (E^) universal [[37], 40 1. In the microscopic 



regime, Leutwyler and Smilga |40| further observed that when ( |45| ) is con- 
verted through Z n (M) = J d6 /2-Ke~ i6n Z(6 , M) to the n-vacuum state, the 
terms of order (VMT,) k in the Taylor expansion of Z n are in one-to-one 
correspondence with inverse moments of powers of the eigenvalues of the 
Euclidean and massless Dirac operator, iIf[A](j)k = \k[A\<p, in the external 
gluon field A. A typical example is the first moment pel 



/ 1 \ V 2 T? 

\?WI/ n = 4(H + jv / ) (46) 

where the averaging is over the quenched gluon field, excluding the zero 
modes (primed sum), in the large volume V limit [| In the n-state, the 
number of Dirac eigenvalues near zero is decreased with increasing wind- 
ings \n\ and flavors Nf, due to the accumulation of zero modes (exclusion 
principle). More importantly, ( |46| ) implies that the spectrum near zero in- 
volves a huge number of eigenvalues ~ 1/U, m comparison to \ j\fV in 
free space. The onset of spontaneous chiral symmetry breaking is followed 
by a large accumulation of eigenvalues in the microscopic region near zero 
virtuality, a point first made by Banks and Casher [43]. 



3.2. Universal Spectral Fluctuations in Vacuum 

Random matrix models share ([45|) with QCD at tree level (macroscopic 
regime) and to leading order (microscopic limit), so they also share these 
moments. This is also true for the effective models announced earlier, such 



as the instanton model [44| and the Nambu-Jona-Lasinio model [pq ]. Unlike 
other models, however, random matrix models provide a simple access to 
a chiral fermionic spectrum where the properties following from the spon- 
taneous symmetry breaking are properly encoded at zero virtuality. (The 
same features can be checked to hold for the Nambu-Jona-Lasinio model or 
the instanton model on an Euclidean hypertorus, provided that the ultra- 
violet contributions in the former are removed.). 

The set of all moments introduced by Leutwyler and Smilga follows read- 
ily from a mean spectral density ^(A) which follows from (filf) as suggested 



8 The ultraviolet singularities of (^) are subleading in V EJ. This justifies the decou- 
pling assumption between soft and hard modes for the inverse moments in the microscopic 
limit. 
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in |40|] . Indeed, for Nf = 1 (see footnote) and in a fixed n- vacuum pC)[ 



1 oo J- 

±V n (\) = Y,6(\V±^) (47) 



k=0 



where £k,n are the positive zeros of J n (£,n,k) = with £ nj £ ~ 7r(/c + n/2 — 
1/4). For V — ► oo and x = AU fixed, (^) reproduces all of the sum 
rules discussed in [40|. This limit is referred to as microscopic, and is best 



taken using random matrix models [ 42 , Q . For continuum four dimensional 



functions [ 33 1 



QCD in the 0-vacuum state, the microscopic distribution is given by Bessel 
xT? 

u s {x) = — [J 2 Nf {T,x) - J Nf+1 {Zx)J Nf ^(Zx)) (48) 



with a node at zero in agreement with (46) and higher moments. The 
larger Nf, the wider the node near zero. Eq. ( f48|) captures the onset of the 
spontaneous breaking of chiral symmetry near zero virtuality as suggested 
by ( f4"7| ) . This result was extended to the n- vacuum state in |45| . 

Although the microscopic distribution was constructed using a Chiral 
Gaussian Unitary Ensemble (ChGUE), it was postulated in [^] to be uni- 
versal, and therefore shared by QCD itself. There, it was also argued that 
for sufficiently random gauge configurations, the distribution of eigenval- 
ues is a priori constrained by the central limit theorem in the microscopic 
regime. This assumption, however, is not needed. Any chiral model is 
bound by chiral power counting to reproduce fl45| ) and hence (|47|). This 
includes chiral random matrix models with an arbitrary polynomial weight 
consistent with unitary invariance (unbroken chiral symmetry) , and positiv- 
ity. Non-Gaussian weights bring about corrections to ( f45f]47| ) that vanish as 
V — ► oo. Hence ( |45| ) uniquely constrains both the diagonal and off-diagonal 
parts of the n-point eigenvalue-correlators in the microscopic limit, forcing 
them to be essentially unique (universal). This point has been explicitly 
proven in [46]. 



In fact chiral symmetry through (^5|) also constrains the fluctuations 
around ( f47D to be generic, since the finite size corrections in 1/V are likely 
to be subleading compared to the bulk distribution, in the regime where 
| (qq) | ~ MVY? . These fluctuations are best analyzed by unfolding the 
spectra ||, that is reexpressing the distribution of states around the mean, 



say (47), and calculating the pertinent spectral statistics. We believe that 
this is what transpires from the detailed calculations in [^7| for noncontin- 
uum QCD. Our observations do not extend to the regime where 1(^)1 ~ S 
except if chiral power counting is only enforced at tree level, or if the volume 
corrections are accidentally small. 
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From our analysis, it is clear that (|45|) is a result that holds for quenched 
and unquenched spectra of the Dirac operator, provided that the associ- 
ated ground state breaks spontaneously chiral symmetry in the macroscopic 
limit. This is the case for QCD, since chiral symmetry is believed to be spon- 
taneously broken even in the quenched case. The effect of the unquenching, 
is to provide for chiral loops etc. , but these effects bring about corrections 
to that vanish exponentially in the large volume limit. This obser- 
vation confirms a detailed analysis of the effects of unquenching on the 
microscopic spectral density [48]. Changes in the microscopic spectral den- 
sity are induced in the transition regime MVT, ~ 1. This is expected from 
chiral power counting since the loop effects are no longer suppressed in this 
regime, e~ ^ M ~ e _1// ^ ~ 1. In this case (f45|) receives contributions from 
loops in the infrared, and is no longer unique. This is the regime where the 
dynamical content of the theory is tested, and where lattice and continuum 
physics are likely to depart. 

We note, that all our arguments apply to continuum QCD. On the lat- 
tice, the issue of chiral symmetry is blurred by the fermion doubling problem 
and the nature of the Lorentz structure on a discrete mesh (staggered formu- 
lations). Nevertheless, some of the above concepts may be extended to this 
case as well. However, we expect various microscopic universalities to set in 
between the strong and weak coupling limits. Indeed, for staggered quarks 
in strong coupling, the exact (but bare) chiral symmetry is U(l) x U(l), and 
is only expected to turn to U(Nf) x U(Nf) in weak coupling. The strong 
to weak coupling limit is followed by a reordering of the quark eigenmodes 
near zero virtuality. This point can be analytically investigated along the 
lines put forward by Gross and Witten p^], and checked numerically using 
present lattice formulations. The true continuum limit is only recovered in 
weak coupling with full Lorentz and chiral symmetry restored. 

Finally, it is worth mentioning at this stage that in odd-dimensional 
spaces the issue of chiral symmetry is more subtle. Nevertheless, the spon- 
taneous breakdown of a global form of chiral symmetry is possible, and 
an analogue effective action to (^5|) for QCD in three dimensions has been 
worked out, together with its microscopic distribution |50t| . This is an im- 
portant point, since the nature of matrix models in QCD encompass the in- 
stanton issue from which they were originally concocted (see below). There 
are no instantons in odd-dimensional spaces. 



3.3. Universal Spectral Fluctuations in Matter 

Under changes in the external parameters such as iVy, N c , M, 8, tem- 
perature, quark chemical potential, QCD may undergo phase changes. If 
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we were to assume that some of these phase changes obey the general lore 
of universality, then at the critical points the transition is dictated by the 
symmetry of the order parameter (again the nature of the coset space) 
and a number of relevant operators as originally discussed by Pisarski and 
Wilczek J5l| ] . While the relevant operators characterize the structure of the 
'effective potential' (radial variable), the nature of the coset in the weak 
field limit dictates the universality of the spectral fluctuations of the QCD 
Dirac operator in matter. For instance, a phase change induced by the 
temperature T would yield around the critical point to {(3 = 1/T) 

Z(9, M; r) ~ / m ^(t^m,*) (4g) 

where £3 involves polynomial operators [] 

£ 3 (T, 6, M, $) = Tr(e ie/N fM$ + ^mV^) + g (T) Tr(&§) 
+ 5 i(T)TV($t$)2 

+ ^(^(det^ + det^) + 0($ 6 ). (50) 

Here $ is a complex Nf x Nf matrix, V3 a dimensionless parameter com- 
mensurate with the three volume, and O are contributions from higher 
dimensional operators as well as non-zero modes. go(T) is negative below 
and zero at the critical temperature. The standard and nonstandard chiral 
random matrix models formulated so far are consistent with ( |50D near the 
critical point 

We observe that near the critical temperature, the relevant operators 
determine the value of \ §\ ~ E(T) at the saddle point in the radial direction. 
In the weak field limit, this value is M independent, and as a result (|4"9]-|50|) 
reduce to © with the substitution E E(T), V /3V 3 , and SU(N f ) 
U(Nf). In the microscopic regime (3V^MT 1 {T) ^ 1, the vacuum sum rules 
are unchanged provided that the concept of an n-state remains valid (the 
9 integration merely bringing an M-independent overall normalization for 
53 7^ 0). This is the case for the chiral random matrix models discussed in 
plf . These arguments are in agreement with detailed calculations in the 
case Nf = 1 and g% = 53]. Following on our arguments in section 3.1.2, 



we would also expect the fluctuations in the infrared part of the bulk of the 
spectrum to be universal, and follow from chiral random matrix models. 
This can also be checked by studying the spectral statistics of the unfolded 
spectra at finite temperature. 

Clearly, our arguments apply to other phase changes provided that the 
analysis is maintained within the spontaneously broken phase (approach 

9 The possibility of phase changes induced by non-polynomial operators, e.g. see [^2[ , 
will not be discussed here. 
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from below), for which the concept of an order parameter such as $ makes 
sense [3- In this sense there is a dichotomy in universality: the universal 
spectral oscillations of the QCD operator in the infrared are tied to the 
universality of phase transitions at the critical point. What is amazing, 
is the fact that the spectral oscillations remain unaltered below the phase 
transition, thanks to the universality of phase changes. 



4. APPLICATIONS of BLUE'S FUNCTIONS 

Having gone over the universal aspects of spectral oscillations in the 
microscopic limit of effective QCD, we now return to some aspects of the 
macroscopic limit, where random matrix models are regarded as effective 
models beyond the contribution of order M and infinite volume. In this 
regime, the 1 /N approximation is an interesting way of organizing the ex- 
pansion (modulo exceptional points), where subtleties may arise from dif- 
ferent limiting procedures (large volume, small masses, quenched, etc.). 

For effective QCD, chiral random matrix models offer a useful alternative 
to conventional power counting as is the case for the U(l) and CP problem, 
and provide a simple framework for discussing phase changes driven by 
conventional and non-conventional universality arguments. Although they 
do not follow the general lore of conventional chiral power counting, they are 
still useful in the infrared for a qualitative description of bulk observables 
for which soft and hard modes decouple (see below). For models of quantum 
mechanics with or without dissipation it is a useful mean for implementing 
mean-field-like approximations. We now proceed to illustrate these points. 



4-1. Chiral Random Matrix Models 

Chiral random matrix models follow naturally from the O-dimensional 
reduction of chiral four-fermi models in even dimensions. Specifically |3l[ 

C = q{iMe ie ' N fq R + q ] R iMe~ w / N f q L + -Lgt^t^ (51) 

where dim qu t L = ^f n +,~ characterizes the dimension of the Grassmanian 
space, and n + + n_ = V . Other 'multi-quark' interactions are possible 
provided that they transform as singlets (1,1) under SU(Nf) x SU(Nf)f^ 

10 For phase transitions in terms of Nf and N c , we should think about interpolating in 
these parameters. Also at finite chemical potential the microscopic limit has to be defined 
with care. 

11 For Nf > 2 with the U(l) broken by determinantal-type interactions. The present 
arguments hold for Nf = 1 as well |pl| . 
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Chirality is the only remnant of the even dimensionality of the original 



space. For odd dimensions, we refer to [50|. In the macroscopic regime, these 



models owe their essential physics to the instanton liquid model (from which 



they were originally inspired) |3^, 41, 42, [54]], and in the microscopic limit 
to the uniqueness of (f45|), following from ( [51] ) by bosonization W = Hlhr- 

Indeed, from present cooled lattice simulations, it appears that the bulk 
features of a number of hadronic expectation values and correlations can be 
well-described by an ensemble of random topological structures (singular 
instantons) where the quark and antiquark zero modes play an important 
(perhaps) dominant role. The hopping between the zero modes in this com- 
plex environment can be assumed to be random with a typical strength set 
by the chiral condensate in the vacuum. This is a non-controllable approx- 
imation, and in general causes most of the results in the macroscopic limit 
to depend on the nature of the random weight. If we were to assume that 
the QCD vacuum is sufficiently random, then there may exist an infrared 
Gaussian fixed point that will cause the weight to truncate to a Gaussian, 
in the macroscopic limit 0- 

Since instanton models are fair effective models 
to the QCD vacuum and its low-lying excitations, we may expect chiral 
random matrix models to be a fair (albeit cruder) effective model for bulk 
observables with weak ultraviolet sensitivity. Their advantages being their 
solubility and clarity in the thermodynamical limit. 

The bosonized version of the four-fermi model ( |5ll) yields W = qnqi, 
that is 



Z{9, M) = W [dW] e -^TrV(wt) . e -^y 

n± 



Here W is a complex random asymmetric n_ x n + matrix. The two diagonal 
blocks have implicit size n=p x n=p. The weight for the random potential is 
assumed Gaussian for simplicity. The widths of the Gaussians: a, x* i K * 
are related to the quark condensate, the topological susceptibility and the 
compressibility of the quenched ensemble, respectively. They have been 
discussed in |56|| . Here 2N = n + + n_, k ± x = 2(n± — (n±)), Nf is the 
number of flavors, and V a dimensionless parameter commensurate with 
the original four volume. For most cases we will assume 2N/V = 1 in the 
thermodynamical limit unless stated otherwise. 

The model (||) has been studied extensively ||, [5], g|, @> USUI 



and this is why we have chosen to analyze it in these lectures. It is minimal, 

12 The existence of an infrared Gaussian fixed point in the macroscopic limit may be 
addressed using the renormalization group method discussed in |pq|. 
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in the sense that it follows from the instanton model to the QCD vacuum in 
the zero-mode sector and infinite wavelength limit |^1J |. To account for the 
presence of near zero modes, along with the topological (exact) zero modes, 
a simple variation of (52) is just 



Z(8,M) = J2 C(n+,n_) f dRe-t Tr(RRt) detQ 

n± ■> 



(53) 



with £(n + , n_) a pertinent measure in the quenched limit, and Q = D + R. 
The random matrix R has rectangular entries 



R 



/ 
At 



A 

o 

T L 

t 



a 



V r R o 



(54) 



/ 



and the sparse and deterministic matrix D has square entries, 

\ 



D 



f iMe ie / N s 

iMe~ ie / N f 

iMe ie / N f 

\ iMe- w / N f J 



(55) 



Each of the four rows and columns in the above matrices (54-55) have dimen- 
sions NfN, NfN, Nfn+, Njri-, respectively. A characterizes the hopping 
between near zero modes, a the hopping between the zero modes, and T 
the cross-hopping between zero and near zero modes. For T = the near 
zero modes decouple from the zero modes in this model. This non-standard 



model is perhaps more in line with the recent QCD studies on the lattice [58|. 



Since (53) share the same (45) with QCD in the vacuum, it fulfills the same 
microscopic sum rules, and hence the same Bessel oscillations |33|]. It is 
amusing to note that a similar class of random matrices (ring matrices) was 
actually discussed by Brezin, Hikami and Zee in some generic models of dis- 
order |59| , where Bessel oscillations identical to the ones following from ( |52| ) 
were derived in the microscopic regime, for the case of symmetric matrices 
with n + = ra_ = N and Nf = 1. 



4.1.1. Decoupling and Spectral Densities 

13 We are restricting the discussion of j3l]] to the vacuum with a 8 angle. The inclusion 
of matter in (^) only affects the near zero modes, ensuring exact zero modes throughout, 
as expected from QCD. 
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We note that both (52) and (^) are chiral in the sense that the argument of 
the determinant anticommutes 75 = diag(l n+ , —l n _), in the massless case. 



For simplicity we return to (52) and restrict for the moment to the case 

om 

(56) 



Nj = 1. For 9 = and fixed k, the spectral properties of ( p2\ ) follow from 

1 \ 




m iW 



Using (||), we have (m — > 0) 

where the density of states v{\) relates to the discontinuity of G(z) through 
(||) with m = 0. For the quark condensate to be non-zero, the average 
number of quark states near zero virtuality has to grow with the size N. 
fl57|) means that the level spacing near zero is of order 1/N in comparison to 
1/yiV in free space. Chiral symmetry brings about huge correlations near 
zero virtuality. The relation (|57| ) was originally established by Banks and 
Casher ||43|. It is also the first of the infinite sum rules discussed recently 



by Leutwyler and Smilga [40]. 



Since the spacing near zero is of order 1/N it follows that in the micro- 
scopic limit iVA ~ 1, the density of states is solely given by the soft modes, 
as the hard modes (essentially free) are separated to infinity. In this limit 
the spectrum near zero is blown up to cover the whole energy scale, causing 
the hard modes to decouple by hand. The random matrix assumptions of 
keeping only the soft modes become legitimate. In the macroscopic limit, 
the hard and soft modes couple, and the decoupling assumption is in general 
not justified. However, for a number of infrared sensitive parameters this 
coupling is weak. For instance, although (|57]) diverges quadratically, the 
divergences are readily seen to be of the form C\m + C2"7- 3 by performing 
a double subtraction. Hence, they do not affect explicitly the value of the 
quark condensate in the chiral limit, as expected for an order parameter. 
They only affect it implicitly through the occurrence of anomalous dimen- 
sions, which is a mild logarithm dependence on the renormalization scale. 
These arguments carry on to matter, thanks to the non-renormalization 
theorems [pCf] . 

The averaging over the rectangular matrices in ((5^) could be performed 



either by using diagrammatic methods [57], or by using the multiplication 
law for Blue's functions by rewriting the rectangular matrices as square 
random matrices multiplied with projectors (see next section). In both 
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cases, G(z) takes the generic form 



G = gl + 5575 (58) 
for fixed n± with 1 = diag(l n+ , l n _). Specifically 

TrG = g + xg 5 = l(z - ^ z 2 - 4 + 4x 2 /z 2 )| 2 ^ m 
Tr G75 = xg + g 5 = --\ z = im (59) 

for the chirality even and odd resolvents with asymmetry x = n_)/2A r . 
The corresponding spectral densities are 

u + (x,x) = M^) + ^\A A2 - A -)( A2 - A +) 

v-(X,x) = x5{\) (60) 

with A| = 2 ± 2Vl - x 2 . For x = we recover Wigner's semi-circle. In 
particular, ^+(0) 7^ in agreement with (p^) and the general lore of spon- 
taneous breaking of chiral symmetry in the infinite volume limit. Note that 
i/_(0) measures directly the difference between the left-handed and right- 
handed quark zero modes. It follows from the rectangularity of the chiral 
random matrices, and is a schematic version of the Atiyah-Singer index 
theorem lEII]. 



4.1.2. Screening of the Topological Charge 



The quenched version of ( |52j ) shows that the topological charge is fixed by 
the width of the Gaussian, that is x ~ X* ~ m ° an d finite in the chiral 
limit. In the presence of quarks, the Gaussian distribution gets correlated 
with the fluctuations in the asymmetry of the matrices. As a result, the 
topological charge % ~ m an d vanishes in the chiral limit. What causes this 
is the occurrence of terms of order x/m ~ 1/Nm in the large N analysis 
which makes the limits N — > 00 and m — > non-commutative. 



To analyze this consider the vacuum partition function (52) for fixed 



asymmetry x = x/2A r , equal quark masses and zero 6. For z = im, we have 

d z log Z x = -2iNN f TrG (61) 



where TrG follows from (|5^). Z x follows through the integration law ( |34| 
or (|35|) for Blue's functions. However, for the screening of the topologi 
cal charge we only need to retain the terms of order x 2 ~ 1/N in G{z) 
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since the Gaussian forces x ~ v N. With this in mind, we have G(z) 
—x 2 /(z 2 Vz 2 — 4), so that integration of Green's (Blue's) function yields 



log Z x = log Z x=0 + — X 2 + 0(x 3 ) • (02.) 



The integration constant was fixed from the asymptotic condition on Z x . 
Inserting Z x into (V = N) 

r -x 2 

Z(0, M) = dx e^iv z x (63) 
yields the following substitution 

Vx* - Vx* + E — r^= l/x * + £ 7i (64) 

i=i mi(J mf + A + mi) i=\ 

after reinstating the flavor dependence. For iVj = (quenched) the topolog- 
ical susceptibility is just given by x* ~ m °> while for Nf ^ it is screened, 
X ~ m <C 1. The fluctuations of order x 2 induced by quark matrices of 
different sizes are crucial for this screening. This mechanism is generic [S2|, 
and holds for more realistic models of the QCD vacuum such as the instan- 
ton vacuum model [ S3 1 . 



4.1.3. U(l) Problem 



The non-vanishing of the unquenched topological charge is usually enough 
to show that QCD solves its U(l) problem. This is also true in the present 



chiral random matrix models. Indeed, at finite 9, (52) allows us to define 
a simple 'Ward-identity' between the topological susceptibility Xi the pseu- 
doscalar susceptibility Xps an d the quark condensate [pTp . For 8 = and 
equal quark masses, this identity is just 

2m m 2 

X = "^2* TrG " J^Xps ( 65 ) 

in total analogy with the QCD Ward identity 

2 c 

iXtop = -j-p (W) + J^2 d 4 x(T*tf;ij 5 i;(x)il)i-f 5 ^{0)} (66) 
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where the pseudoscalar susceptibility is defined by the integral in ( pq ) . The 
resolution of the U(l) problem in QCD stems from the fact that the un- 
quenched topological susceptibility xtop ^ 0, for otherwise (|66|) implies that 
the powers of m matches only if the pseudoscalar susceptibility develops an 
infrared sensitivity in the form of 1/m, much like the pion susceptibility. In 



the chiral random matrix model, the analogue of (6€) is 



X P s = --^(trQ-SQ" 1 ^} + -^(trQ-VstrQ-Ss)^. (67) 



where Q denotes the argument of the determinant in (52). From the pre- 
ceding discussion, it follows that 

Nf 1 1 N f 

which is finite in the chiral limit. For illustration, we show in Fig. ^ the 
behavior of the various terms in (|68|) following from ensemble averaging 
large samples of random matrices of different sizes N versus the analytical 
results (solid curve) in the large iV limit. 

The remarkable cancelation of the the 1/m 2 terms in ( |68] ) between the 
trace-term and the trace-trace term, implies that the disconnected term is 
as important as the connected one in this channel, contrary to conventional 
wisdom. The reason is the strong infrared sensitivity in this channel as 
emphasized above. Crucial to our argument is the issue of the large N limit 
or the thermodynamical limit since V ~ N in our analysis. The above issues 
are important, and should be retained in a lattice analysis of full QCD, as 
they teach us valuable lessons on the subtle interplay between fluctuations 
in finite size lattices and the onset of screening. 



4.1.4. 9 Vacua 



For 9^0 the determinant in (52) is complex. Simulations in QCD with 
complex fermionic determinants are still elusive, owing to the occurrence 
of random phases. In simplified models, we refer to the recent numerical 
studies in |j4|. For this problem, chiral random matrix models may be 
effective in addressing some of the underlying issues, especially those related 
to the thermodynamical limit. Having said this, we note that ( |52"| ) is periodic 
of period 2tt in 9, provided that the sums n± span over n- vacua with n even 
and odd. We also note, that in general the matrices sampled in ( |52"D are finite 
in size, thereby making the issue of the 9 dependence of the vacuum state 
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Figure 2: Topological (left), pseudoscalar (middle) susceptibilities and 
anomalous Ward identity (right) for three equal flavors. The numerical 
simulations were carried for for fixed 2N = n+ + n_ = 60, 120, 200. Solid 
lines are analytical results. 



subtle in the thermodynamical limit pq ]. In a nutshell: the precision of the 
calculation and the thermodynamical limit are intertwined, and erroneous 
results may emerge if due care is not taken |6f| . The large TV" analysis we will 
present below has been found to be consistent with the above observations 

With this in mind, and to investigate the behavior of the vacuum energy 
versus 8 we now proceed to estimate lnZ(#, M)/N . If we were to assume 
that the range of resummation over x = n + ~ n ~ is infinite, for N = n++n_ 
peaked with a small but nonzero spread, then the Gaussian integration in 



( f52|) can be readily performed. Standard bozonization 6J] gives 



hxZ(6,M) = iVTrln U±±J^j + ^ (rrln^^ (69) 

where z = diag mj e id / N f, and P is a complex Nf x Nf matrix expressed in 
terms of the TBy's and 9 through the saddle point equations (see below). In 
the vacuum, P = diag pje^ s ^ N f~^' , so that the unsubtracted free energy 
associated to (|69|) reads 

F(6, M)= ( ln Pj ~ ln (P 2 j + ZrrijPjCOS 4>j + mf 




(70) 

From ([7(]) it follows immediately that if the quarks masses are large, 
that is nij ^> pj ~ 1 at the saddle points, then the 6 dependence of the 
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unsubtr acted free energy is just that of the quenched limit, that is F ~ 
X*6 2 /2. For small masses rru <C x* ~ Pj ~ 1 (the scale of spontaneous 
symmetry breaking), then ([70|) simplifies to 







F(e,M)= Y,[p 2 3 


— Inpj - 


j=i \ 




2N \ 


( Nf 

*-E 

V 3=1 



rrin 



Pj 



E 



-sine 



; Pi 



+ 0(m 2 ). (71) 



The saddle point solution in the p's decouples and gives pj = 1 + 0(m) 
The saddle points in the </>'s give 



3=1 



mismcpi 



(72) 



in agreement with results derived by Witten [39] using large N c (number 
of colors) arguments. That they are reproduced by ( |52| ) is reassuring. This 
also means that they should also be satisfied in the instanton model as 
well as some variants of the Nambu-Jona-Lasinio model with proper U(l) 



breaking [66]. 

For Nf = 2 these equations can be readily solved to give 



sin0i 2 



m2 isin# 



777^ + 771-2 + 2777-17772COS 6 



(73) 



As a result Z(6, M) is a smooth 2tt periodic function except at 8 = tt and 
mi = 1^2 where both the numerator and denominator vanishes in (|73|). In 
this case, the free energy is again 2tt periodic but with a cusp at 6 = tt. 
The cusp reflects on a first order transition caused by the coexistence of two 
CP violating solutions. As a result CP is spontaneously broken at 9 = tt in 
agreement with earlier arguments [[3^, 67 , |68| , |69| ] . 

For Nf = 3, the explicit solutions to (72j) are in general involved analyt- 
ically. However, at 9 = tt the analysis simplifies. There are two classes of 
solutions. In general, a unique solution that is CP symmetric, and a doubly 
degenerate solution for 77717772 > 7773] m\ — 7772] for which CP is spontaneously 



broken, again in agreement with a number of analyses [39, 67, 69] 



4.1.5. Finite Temperature 
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The formalism of Blue's functions is very well suited to investigate the 
random matrix models of the QCD Dirac operator for finite temperatures 
and/or chemical potential. Indeed, the modifications in the medium corre- 
sponds to replacement of the argument of the determinant (modulo trivial 
factor of i) ( |52] ) by 

where [i is the chemical potential for "quarks" and fi = uj n l n ® ljy. Here 
uj n = (2n + X)tvT are all fermionic Matsubara frequencies. The form of the 
replacement in (74) is consistent with the effective models of QCD and the 
Dirac structure [31] near the critical temperature. For simplicity, we restrict 
ourselves to the lowest pair of Matsubara frequencies ±7rT and the chiral 
limit. This corresponds to the model proposed in [f70f |. 

To investigate chiral symmetry breaking in this model we should calcu- 
late the Green's function and, through the Banks-Casher relation, analyze 
the chiral condensate near the critical point. Therefore, we have to "add" 
the deterministic (D) and random (R) "Hamiltonians" , 

ttT \ ( W \ , . 

irT ) + {w^ ) ■ (75j 

In the present form, the Blue's function technique does not incorporate the 
chiral block structure of the matrices, but for the deterministic matrices 



considered here it works correctly [71]. The Blue's function for (|75| ) satisfies 



B(z)-z = B D (z) (76) 

where Bp is the Blue's function of the deterministic piece. By definition we 
have 

z-G(z) = B D (G(z)). (77) 

Now if we evaluate the Green's function of the deterministic piece on both 
sides of this equality we will get Pastur equation |t|] 

G D (z - G{z)) = G D (B D {G{z))) = G{z) . (78) 

The explicit form of the deterministic Green's function is Gd{z) = z/{z 2 — 
tt 2 T 2 ), since the deterministic eigenvalues are ±7rT, so the final result 
reads |^ 

G 3 - 2zG 2 + (z 2 - tt 2 T 2 + l)G-z = 0. (79) 

14 This cubic equation is shared by chiral and non-chiral random matrix models alike . 
The chiral structure shows up only in 1/N. 
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Figure 3: Spectral density u(X) as a function of the deterministic parameter 
e = irT. 



The 'critical temperature' for chiral symmetry restoration follows from the 
the condition B'{G) = 0. Using (|79|) this locates the end-points of the quark 
spectrum at 



1 (4tt 2 T 2 - 1 + V8tt 2 T 2 + 1)§ 
V8ttT V8vr 2 T 2 + 1 - 1 



_ i (4,^-1- Vi^Ti), 

a/87tT V8vr 2 T 2 + 1 + 1 

The condition A2 = corresponds to the situation when the two-arc support 
of the spectrum [— A\, —A2] U [A2, Ai] evolves into a one-arc at [— Ai, +Ai\. 
These structural changes in the spectrum of quark eigenvalues are shown 
in Fig. ||. The critical temperature is T* = 1/V, in units where the width 
of the random Gaussian distribution is set to 1. The break up near zero 
virtuality (A ~ 0) may be generic of a phase change in the system, provided 
that the decoupling between soft (S) and hard (H) modes take place in the 
Dirac spectrum 1 31 ] . This is the case for order parameters as we discussed 



above. 

The set of all critical exponents could be easily inferred from solving 
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the algebraic equation ( |79| ) and calculating the free energy of the system by 
integrating Blue's functions. The result is 

F = G 2 + log^^. (81) 

All critical exponents are of the mean- field- type in agreement with conven- 
tional universality arguments It is not difficult to construct random 
matrix models with non-mean field critical exponents. Using the results in 
| p If one can easily see that the effective potential at the critical point is 
all what matters, much like the effective potential (^5|) at zero temperature 
as discussed in section 3.1. A non-mean field effective potential can be un- 
raveled in the form of a chiral random matrix model by using the inverse 
procedure. An example of a model with first order transition was given in 
©■ 

How realistic are all these statements? This can be seen by comparing 
some of the present results to lattice simulations near the critical point. 



Recent simulations by the Columbia group [73|, have revealed interesting 
insights on the distribution of eigenvalues using the semi-quenched conden- 
sate 

<«>= 2 ™<f rfA A^ < 82 > 

where the distribution of eigenvalues is for Nf = 2 with sea quark masses 
m s 7^ nv. Although at finite there is a sensitivity of (^) to the hard 
modes (see above), the latter is suppressed for small m^. The behavior of 
( |82| ) for staggered fermions on a 16 3 x 4 lattice is shown in Fig. |] (left) versus 
the valence quark mass for equal but fixed sea quark masses m s a = 0.01 
(6 MeV). In the chiral random matrix model, the eigenvalue distribution 
v{\) follows from @ through the substitution ir 2 T 2 -> ir 2 T 2 + m 2 @. 
The result for (|8^) is shown in Fig. [| (right) for m s = 0.1 (10 MeV). We 



have identified the lattice (3 with T through : ir(T — T c ) = {(3 — (3 C ) and set 
Pc = 5.275. There is overall qualitative agreement over several decades of 
the valence quark mass m^. 



4.1.6. Finite Chemical Potential 



The case of finite chemical potential can be analyzed in a similar way, using 
instead 

in \ , ( 



J + U o ) • ^ 
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Figure 4: Semi-quenched condensate versus the valence quark mass for 
two- flavor QCD |73[ ] with a sea mass m s a = 0.01 (left). The result on the 
right is from a chiral random matrix model. See text. 



The deterministic part is now nonhermitian and we have two distinct regions 
(at least in the quenched case) on the complex plane. The Green's function 
in the holomorphic region outside the blob of eigenvalues is given by ( f79| ) 



with the substitution tt 2 T 2 



2. 



G 



2zG 2 + (z 2 + fj, 2 + 1)G 



0. 



(84) 



In order to find the boundary of the nonholomorphic region we may exploit 
the method of conformal mapping via the Blue's functions, and transform 
the cuts of the T ^ case into the boundary (see Fig. |5|) by the transfor- 
mation 

z — ► i w = z - 2G(z) (85) 

where G(z) is the appropriate branch of the cubic equation ( f79|) with the 
formal identification tt 2 T 2 = +/x 2 . The extra factor of i corresponds to the 
rotated Hamiltonian (here the deterministic as opposed to the random part 
is antihermitian) . 

The result of this mapping is in total agreement with a result derived 



originally by Stephanov [75] using a variant of the replica technique for the 
identical problem. The domains encircled by the mapping flS5|) defines the 
support of the complex eigenvalues of (|83|). In this region the resolvent 
G(z) is non-analytic. The analytical resolvent following from (|84] ) defines 
the correct branch in the holomorphic (shaded region), but fails on the 
boundary where the 1/N expansion breaks down. The failure is caused by 
the accumulation of eigenvalues in the plane, and signals large fluctuations 
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Figure 5: Conformal mapping from temperature to chemical potential. The 
one- and two-arc cuts (left) are mapped onto islands (right). The shaded 
region is the holomorphic support of the mapping. 



in the phase of ( |83| ) as well as its correlations [76]. Similar effects are 
observed in quenched lattice calculations |77J . 

This effect is an artifact of the quenching (suppression of the fermion 
determinant) in the calculation of the resolvent which is a singular operator 
This may be overcome by unquenching at the price of much larger 



precision in the numerics, due to the presence of phase oscillations [78], or 
analyzing non-singular operators such as the partition function, which are 
self-quenched in the large N limit [18, ^ll [J4j|. High precision calculations 
for the unquenched partition function Z^(z = 0, m) with Nf = 1 are shown 
in Fig. |6] (dotted line). The solid lines refer to the two holomorphic so- 
lutions of (84). The long arrows indicate the points where the islands of 
non-analyticity in the quenched resolvent cross the real axis. The short 
arrow is the position of the cusp resulting from the intersection of the two 
holomorphic solutions to ( [84|). 

The phase structure of (|33D directly follows from the analytical prop- 
erties of the partition function Zn(z,[i) for general complex z = im with 
Nj ^ and in the 1/N approximation as also noticed in fT"5|] . In our case, 
this can be readily done using the integration law for Blue's functions and 
the holomorphic solutions to (|84"|) [18]. The outcome is simple: Zn develops 
cusps in the z- or ^i-plane that reflect on the two-branches of (84), with 
end-points conditioned by B'iG) = 0, in agreement with the results in ||75|. 



The curves of zeroes obtained in this way are also in agreement with the 
(high-precision) numerics of [79]. The location of all zeroes follow from the 
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Figure 6: Unquenched free energy for N f = 1 at /i 2 = 0.8 (left) and /i 2 = 1.2 
(right). The dots are high precision numerical results, while the solid lines 
follow from the solutions to fl84|). See text. 



analytical condition (|33"D, here 



Re 



G 2 (z)-G 2 (z)+log. 



_^) 

G,(z) + 8 G,(z) 



(86) 



where the labels 



refer to the different branches of the solutions 



The density of zeroes follows accordingly from (37). 

Using similar techniques one could also calculate the effect of fluctuations 
(l/N corrections to the one-point Green's function), with the result [18] 



Z N (z,n) 



o NE 



2 [(D 2 + M 2 ) 2 



{z-Gf]\ 



(87) 



Both Eq and D depends functionally on G to this order. Hence, the infor- 
mation carried by the Blue's functions condition the analytical behavior of 
fl87P. In particular, the location of the zeros of the curly bracket coincide 
with the points at which the 1/N expansion breaks down. These are ex- 
ceptional points that may signal the onset of new scaling regimes for a new 
class of microscopic universalities in the context of non-hermitian random 



matrix models [18 



4.1.7. Phase- Diagram 



For general finite T and finite [i the discussion about the thermodynamics of 
chiral random matrix models is very much model dependent, much like the 
thermodynamics in the Nambu-Jona-Lasinio model or the instanton liquid 
model. It is simply that of 'constituent quarks', with one quark and one anti- 
quark state (two-level) [j3l| , j34f . In the absence of matter, chiral symmetry is 
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Figure 7: Sample distribution of eigenvalues for different fx 2 . The solid lines 
follow from the conformal mapping. 



spontaneously broken by filling the 'negative-energy' state. In the presence 
of matter, the symmetry is restored by depleting the 'negative-energy' state 
and filling the 'positive-energy' state with equal weight (finite T), or when 
the quark chemical potential reaches the 'positive-energy' state (finite fj,). 
In the first case the transition is second order, while in the first case it is first 
order. Diquarks effects in the present chiral random matrix models are down 
by 1/N, and hence they decouple in the large N limit. They are inessential 
for the thermodynamics in these models. Although unwelcome for N c > 3 
QCD because of color confinement (zero triality), they are dominant for 
N c = 2 QCD where they play the role massless baryons. 

Having said this, we now consider (|74|) with infinitely many Matsubara 
modes. This is essential when departing from the critical points, or in 
taking the temperature to zero for fixed \x [|31f| . In the chiral limit, the quark 
condensate \ {qq)\ = P* (spectral density ^(0)), follows from the saddle point 
equation to the partition function associated to (74), that is 34] 



= 1 ' ' 



I + e (|P.+m|-Ai)/T i + e (|P*+m|+ M )/T • 

The solution to ( |88|) for m = is shown in Fig. ||. For fi > 1 /6S In (2 + \/3), 
the transition is first order and disappears at //* = 1/2S. At fi = 0, the tran- 
sition is second order (mean-field) and sets in at T* = 1/4E. The mismatch 
with the = 1/tt quoted in section 3.1.6 is due to: a different normalization 
for S and an infinite number of Matsubara modes. The behavior shown in 
Fig. is overall consistent with results in four dimensions using constituent 



quark models, e.g. [35]. We refer to pT| , 34 for further details on some of 
these and other issues in relation to phase changes. We note that the result 
for n = was also reproduced by [^fl] using a similar random matrix model. 



15 Alternative chiral random matrix models with diquark states to order N° are straight- 
forward to write down. They will be discussed elsewhere. 
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Figure 8: \{qq)\ from (Eq) with X = 1/4 and m = 0. 



The effects of the inclusion of more Matsubara modes was also investigated 
in §3|. 



4-2. Strongly Nonhermitian Ensembles 
The theory of open quantum systems plays a key role in many areas 



of physics and chemistry [83]. Microscopic treatment of dissipative evolu- 
tion involves hermitian Hamiltonians. However, the use of non-hermitian 
Hamiltonians has been found to provide a convenient way to give a reduced 
description of the system. In this case, the hermitian part of the Hamil- 
tonian refers to the free undamped dynamics whereas the nonhermitian 
remainder describes the damping imposed on the system by some external 
source. In this section, we will show a few applications of the Blue's func- 
tions for dissipative systems, where the effective Hamiltonian is obtained by 
the standard Wigner-Weisskopf [84] reduction by partition. Reduction by 



partition follows from dividing the whole Hilbert space into two subspaces, 
and then the Hamiltonian is integrated over one of them. In this way one 
reduces the eigenvalue problem of high dimension to a lower one, and the 
resulting nonhermicity of the effective Hamiltonian comes from the "leakage 
of probability" due to the truncated vectors of the original Hilbert space. 



4.2.1. Open Chaotic Scattering 
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We start from non hermitian random matrix model introduced by Mahaux 
and Weidenmiiller |85|] for resonance scattering. In brief, a quantum system 
composed of (N — M) closed and M open channels can be described by the 
effective scattering Hamiltonian 

H = H R + H T = H R -igT , V = AA T (89) 

which is N x N dimensional. A is an asymmetric N x M random matrix, 
Hr is random Gaussian Hamiltonian and g an overall coupling. Unitarity 
enforces the form of T used in ( |89[ ) . The matrix elements A^, characterize the 
transition between the (N — M)-internal channels and M-external channels, 
and are assumed to be independent of the scattering energy. For M = 0, the 
Hamiltonian is real and the spectrum is bound. For M^O, the Hamiltonian 
is complex, with all states acquiring a width. In what follows, we will solve 
the system by "decomposition" into subsystems and consecutively apply the 
methods presented above for Blue's functions. 

First consider the case of an NxN real symmetric product matrix = 
AfAf. Note that the Green's function (for even random potentials) for the 
square of a square matrix is related to the Green's function for the unsquared 
one via G(z 2 = w) = G(z)/z, due to the identity 

From (|9C| ) and (||), we infer the resolvent for the square of the Gaussian 

G* = i(l- f^) . (91) 

Second, let us observe, that the case of rectangular N x M matrices 
with N, M — > 00 but m = M/N fixed, follows by truncation using the 
projector p6| 

P = diag(l 1 _ ; _ I l, (V^) (92) 

M N-M 

that is T = T S P. We recognize the problem of "multiplication" of the 
random matrix V s by the deterministic projector P, and therefore could use 
the multiplication law for the Blue's functions, hence the S-transform. The 
construction of the resolvent for the projector is straightforward 

= ^ Tr 7^p = m YZ~i + C 1 - ■ ( 93 ) 
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Using the relations for "multiplication law" from the previous part we obtain 
the S-transform 



Sp(z) 



1 + z 

m + z 



(94) 



Similarly, the S-transform for the Green's function for the square of the 
Gaussian (|9l|) is 



S A 2 



1 



l + z 



The product reads 



S = Sjg • Sp 



1 



m + z 

so inverting the order of reasoning we get for the resolvent 



G(z) 



1 — m 1 

2z + 2 



1± 



1 — m 



z 



1 



Am 



(95) 



(96) 



(97) 



Note that the 1/z terms (first term and a term under the square root) 
represent (JV — M) zero modes of the matrix T. 

Now we could "add" the Hamiltonians Hp and —igT, by constructing the 
generalized Blue's functions. Let us first write down the hermitian Blue's 
functions for Hp and The first Blue's function, coming from inverting 
(|) is 



B 1 (z) = z + -. 

z 



The second one is the functional inverse of (|97|), therefore 

mg 



Bo z) 



1 

+ -• 



1 - gz z 

Then the generalized Blue's functions read 

Bi(Z) = 2 + -^ 

B 2 (Z) = m(l-^Z)- 1 g+^ 



(98) 



(99) 



(100) 



where the coupling g is a matrix g = di&g(—ig,ig). The signs in the ma- 
trix coupling reflects the fact, that the complex conjugate "copy" interacts 
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Figure 9: Comparison of the numerically generated distribution of eigen- 
values for few sample couplings in the model (|89| ) with the analytical re- 
sults. The solid lines represent the boundaries following from the condition 



G gq - = 0. 



with opposite coupling constant. Then, the generalized addition law for the 
Blue's functions leads to the final matrix equation for the Green's function 

Z = m(l-ggy 1 g + g + ^. (101) 



Nontrivial (nonholomorphic) solution to the matrix equation (101) provides 
full information on the eigenvalue distribution. The condition G qq - = 
defines the boundaries of the eigenvalue distribution, and the divergence 
of G qq (Gauss law) gives the desired spectral density of eigenvalues on the 
islands. Instead of presenting the explicit formulae, we show the figures 
with few sample solutions (figure ||) . 

The results obtained above were obtained also by other methods []87| , so 
we do not need to dwell on their physical content. However, we encourage 
the reader to compare the present derivation with the derivations of the same 
result based on the replica or supersymmetric methods (see e.g. |88], |89fl)- 



4.2.2. Bridging via Dissipation 



We conclude this brief overview by another example, where analytical results 
were obtained by direct application of the generalized Blue's functions. This 
is the model of an effective two-level system, coupled to a noise reservoir 



39 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 



XXX 

Figure 10: Comparison of the numerically generated distribution of eigenval- 
ues for few sample parameters of the model ( |102| ). The solid lines represent 
the boundaries following from the analytical condition G q q = 0. 



of the random type [90|. Such systems could e.g. serve as model systems 
for electron or proton transfers in condensed phase reactions in the field 
of chemical physics. The deterministic part of the Hamiltonian is given by 
an effective two-level system with N electrons, building up two close energy 
states with half of the energies equal to e and other half equal to — e. For the 
dissipative part, we take H-p from (|89|), therefore the effective Hamiltonian 
reads 



H = H D - igAA ] . 
The deterministic Green's function for two levels is 

c i<J_ + T i_) 



2z 



z + e' 



(102) 



(103) 



therefore the relevant generalized Blue's function is given by the matrix 
equation 



2 = 1 

2 



(B D {Z) - e)- 1 + (B D {Z) + e 



(104) 



The generalized Blue's function for the noise term is just B2 in (|10C| ). The 
generalized addition law for the Blue's functions leads at once 



-5' 



[Z - m(l - gQ)- l g - e)" 1 + (Z - m(l - gG)~ l g + e)" 1 ! . (105) 



As in the previous example, one easily infers the spectral density of the 
system ( |102j ) from the nonholomorphic solution to the equation. Again, we 
present a figure instead of the explicit formulae (see figure |l(]) . 

We would like to note two "phase transitions" visible easily from the 
change of character of the envelopes. Here the envelope is represented by 
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Figure 11: Eigenvalue distribution in the Hatano-Nelson model for large 
(left), intermediate (middle) and small (right) h. The two dots refer to the 
location of the forking (A) and the maximal imaginary eigenvalue (B). 



the solid curve, obeying the analytical solution to G q q = 0, here given by 
the 8th order polynomial in x and y ( with z = x + iy). The first phase 
transition (figure |To|(left) — » ^(middle)) corresponds to the emergence of 
the noise-induced "bridging states". The second phase transition, present 
at large coupling constant (figure |i~0|(m iddle) — > |i~0|(right)) corresponds to 
the "collectivization of the spectra" |9lf|. The upper island represents di- 
rect resonances, the lower island represents the (long-lived) coherent state, 
caused by the explicit presence of zero modes in the dissipative part of the 
Hamiltonian. 



4-3. Weakly Nonhermitian Ensembles 



As a final application of the Blue's function techniques we will consider 
the model of Hatano and Nelson p^] (HN) for nonhermitian localization, 
which has recently attracted much attention ||, |93|, [)4|]. This model is a 
non-hermitian version of the Anderson model, hence quantum mechanical 
in nature. Nevertheless, the method of Blue's functions will turn out to 
be very appropriate for its treatment, since it is equivalent to the Coherent 
Phase Approximation (CPA) widely used in the treatment of hermitian 
models with site disorder |9F 



In their original analysis, Hatano and Nelson [92| have shown that the 
depinning of the flux lines from columnar defects in superconductors in 
d + 1-dimensions may be mapped onto the world-lines of bosons on a d- 
dimensional hyperlattice. The localization by defects was approximated by 
an on-site randomness on the hyperlattice, and the depinning induced by the 
transverse magnetic field was found to cause a directed hopping of the bosons 
on the hyperlattice, resulting in a nonhermitian quantum Hamiltonian. 
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The HN model (for d = 1) is given by the following Hamiltonian 



N . 

H = H + V A =J2n ( ehc A+i c A + e- h c j A _ lCA ) + V A c\c A (106) 
A=l 

where the c A is a boson creation operator at site A. The lattice spacing is set 
to 1. The hopping strength is te ±h , where h is the strength of the "transverse 
magnetic field" |)^] in units of the flux quantum. Diagonal randomness 
is represented by random variables Va, originally taken from a uniform 
distribution in the interval [—A, A]. For h ^ 0, the eigenvalues are complex. 
As stated earlier, this is a nonhermitian analogue of the Anderson model in 
d = 1 dimension. The eigenstates of this Hamiltonian with an eigenvalue A 

(2) 

such that ImA 7^ are delocalized. For sufficiently large h > h c all states 
are delocalized and form a closed curve that looks like a deformed ellipse (see 

figure |ll](left)). When h is reduced below the eigenvalues at the edges 
snap to the real axis and get localized (see figure |Tl|(middle)). Reducing h 

further below causes the eigenvalues to fill in the whole interval on the 



real axis (see figure 11 (right)). Our aim here is to understand these effects 
quantitatively, using the method of Blue's functions. 

The relevant quantity to calculate is the Green's function 

G ' 2 ' = J? ( tr 7~~h) ' < 107 > 

Its holomorphic part is related to the derivative of det(z — H) outside (large 



z). A close look at (106) shows that before averaging over Va, 



/t\ N r n 
det{z-H) ~ f-J -[e Nh + e- Nh \+ f(z,N,V A ) (108) 

where the function / does not depend on h and has the leading order be- 
havior 

f(z, N) = e Nfouti - zyA) ■ (1 + corrections) (109) 

where the corrections are polynomial. For large N there are two regimes: 
h + log(t/2) > f out or < f out . In the first regime the determinant is 
z-independent, hence G(z) = 0, and we are inside the curves shown in 
|TT|(left) , (middle) . In the second regime the determinant is /i-independent, 
hence G{z) is also h- independent, and we are outside the curves. 

(1 2) 

Once we have f ou t we may calculate the critical parameters h c , for 
which we have full characterization of the eigenvalue spectra. Note that f ou t 
depends explicitly on V A . It may happen that for certain configurations of 
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site randomness Va the / term wins while for others the l h' term domi- 
nates. This leads to a broadening of the curve of eigenvalues. This effect is 
difficult to control analytically, however numerical simulations suggest that 
this broadening is of order 1/yN ||. This behavior is reminiscent of the 
case of weak nonhermiticity discussed recently in J96| |. 



4.3.1. Coherent Phase Approximation 

A direct calculation of G(z) is involved. If we were to make a diagrammatic 
expansion we would find that a lot of nonplanar graphs would contribute. 
Furthermore one cannot do a perturbative expansion around h small or 
big, neglecting the 'hopping parameters' above or below the diagonal of the 
Hamiltonian. This is not surprising, the HN model is just the nonhermitian 
version of the Anderson model in d-dimensions, for which there is no exact 
solution in large N except for the special case of Cauchy randomness (Lloyd 
model) H. 

We have shown in H that by assuming the site randomness and the 
'hopping' part of the Hamiltonian to be free, we can just apply the method 
of Blue's functions to the present problem. In this approximation, we resum 
at once the planar diagrams and all non-planar diagrams with single-site 
rescattering (some sort of a mean-field analysis) . This approximation bears 
much in common with the CPA method as we now show. Specifically, let 

H = H D + H R (110) 

where Ho is the deterministic hopping term and Hr is the diagonal random 
matrix with random entries Va- Let Ge> denote the deterministic Green's 
function: 

G D (z)=tr—^— . (Ill) 
z - H D 

In the CPA the Green's function of the system is given by 

G(z) = G D {z-Z) (112) 
where £ satisfies the condition of 

V A -E 



i „ v , r , ( w 0. (113) 
1 - (Va - Y,)G(z) 

Multiplying by G(z) and adding and subtracting 1 in the numerator we get 

1 



1 - (V A - E)G(z) 



1 (114) 
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then multiplying by G{z) once again, we obtain 



G(z) = G R (~ + e) . (115) 

As a final step let us evaluate the Blue's function for the random part Br 
on both sides of this equation 

B R (G(z)) = -^- + X. (116) 
But this is a direct consequence of the addition formula for Blue's functions 

B(z)=B R (z) + B D (z)--. (117) 

z 

Namely after reexpressing Br in terms of B and Bp we obtain 

3 -B D (G( 2 )) + ^ = ^ + E. (118) 



But now we may use the definition of £ ( |112| ) and get 

z-B D (G D (z-Z)) = X (119) 

which finishes the proof. 

We note that the above arguments do not rely on whether the start- 
ing Hamiltonian is hermitian or not. In the hermitian case our conclusion 



is in agreement with the one of Neu and Speicher |11], who used different 
arguments. Since the Green's functions for the nonhermitian and hermi- 
tian ensembles coincide in the "outer" region, the addition of the Blue's 
functions together with the balance of the free energies as detailed in the 
preceding sections, provide for a concise approximation scheme for nonher- 
mitian localization in the large N limit. We now proceed to implement 
this. 



4.3.2. Nonhermitian Localization 



The Green's function for the deterministic part is just 

r M = J_ v = I l/Vz 2 -f 2 "outside" 

N ^ l z-tcos(2Trn/N + ih) \ "inside" 



(120) 
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where the last equality holds in the large iV limit. This shows that we may 
define the Blue's functions only outside the curve of eigenvalues. There it 
is given by 

B D {z) = sj^+t 2 - (121) 

The Blue's function for the site randomness Va is fixed by the type of 
distribution we select. For the (solvable) case of a semicircular distribution, 
the Blue's function is 

B s g micircular {z) = z + - z (122) 

while for a uniform distribution of width A it is 

B u ™ farm (z) = AcothA* (123) 
and for a Cauchy distribution^ 



Therefore the Blue's function for the HN Hamiltonian is 

B(z)=B D (z)+B R (z)-~. (125) 

So the Green's function (in CPA) would satisfy the equation B(G(z)) = z. 
In particular, for a uniform distribution the Green's function satisfies 



-4 + t 2 + AcothAG - l/G (126) 



while for a Cauchy distribution it satisfies 



J_ +t 2_ il = z . ( 127 ) 

Integrating the Blue's function for the HN Hamiltonian gives the free energy 
f ou t (in this case a direct integration of G would be infeasible) 



f out (z) = AGcoth AG - 1 + log 1 + - log (128) 

for uniform randomness and 



f out (z) = log(* 7 + z + t/(i 7 + zf - V) - log 2 (129) 



16 For ensembles with unbounded moments, the proof of the addition formulae is more 
subtle, cf. H. 
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0.5 1 1.5 2 0.5 1 1.5 2 
h h 



Figure 12: Dependence of A (forking point) and B (top point on the imag- 
inary axis upon lift-off at z = 0) versus h for A = 1 (uniform randomness). 
The solid lines follow from the analytical results using the addition law 
(CPA). 



for Cauchy randomness. 

We will now calculate h^P — the value of h where the curve of eigenval- 
ues starts to appear (i.e. when point B in figure 11 leaves the origin along 
the imaginary axis). This is determined by the condition that the boundary 
between the two phases passes through z = 



/i + log^ =/ OHt (0). 



(130) 



For a uniform distribution (128) the solution of the analytical condition 



is only possible numerically. In figure 12 we compare this solution to 
the exact numerical values for the HN Hamiltonian. The agreement is fairly 
good. For a Cauchy distribution, the critical condition (|13C|) is simply 



sinh hffl 



2 
t 



The same condition was obtained by Brezin and Zee 
and Khoruzenko [93], using different arguments. 



(131) 

], and Goldsheid 
We note, in agreement 



with these authors, that in the case of Cauchy randomness the results are 
exact, much like the Lloyd model in the hermitian case. 



5. CONCLUSIONS 

Voiculescu's free random variable approach has a simple realization in 
terms of Blue's functions, as originally suggested by Zee||. This applies for 
both the R- and S-transforms. In the first part of these lectures, we have 
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gone over the essentials of Blue's function calculus, emphasizing the generic 
aspects which are important in physical settings. The addition and multi- 
plication of free random Hamiltonians translates simply to the addition and 
multiplication of the pertinent Blue's functions. More importantly, their 
extrema are in one-to-one correspondence with the end-points of pertinent 
eigenvalue distributions, and their analytical structure conditions the sin- 
gularity structure of appropriate partition functions. Moreover Blue's func- 
tions provide a surprising method of finding the boundary of the support of 
eigenvalues for nonhermitian ensembles through conformal mappings. 

The free random treatment of matrix models amenable to a 1/N expan- 
sion is usually exact to leading order in 1/N, hence in agreement with the 
treatment based on more conventional methods. The advantage of the for- 
mer being its simplicity. For disordered models which are not immediately 
amenable to a 1/N expansion, the free random treatment provides for an 
approximation that is identical to the CPA approximation. The latter has 
been successfully applied to the Anderson model Q, and its variants [11, 12]. 
The free random variable approach goes beyond Gaussian measures, and ap- 
plies equally well to random matrix models as well as quantum models with 
diagonal disorder ||. 

In the second part of these lectures we have discussed effective QCD 
in light of the spontaneous breaking of chiral symmetry in the vacuum. 
Using power counting both in the microscopic and macroscopic limit, we 
have argued that the partition function in the large volume limit is solely 
conditioned by the mode of symmetry breaking on the pertinent coset span 
by the little group [37, 4C[. The stability of the partition function against 
soft and hard quantum corrections, conditions the universal behavior of the 
spectral oscillation of the QCD Dirac operator both at zero virtuality and in 
the bulk of the infrared region. For chiral random matrix models this means 
simply that arbitrary polynomial weights that are chirally invariant yield 
by power counting the same universal spectral oscillations in the vacuum. 

The universal behavior of spectral oscillations in effective QCD occurs 
in the microscopic limit with light quark masses (for chiral symmetry to 
be relevant) and large enough volumes (for nonperturbative physics to set 
in). This is the limit, however, where the vacuum has no spontaneous 
alignment yet (the chiral condensate is driven by the light quark masses). 
The universality arguments do not carry any dynamical content. Rather, 
they reflect on the kinematical correlations induced by the compactness of 
the coset space associated to the would-be Goldstone modes, in the presence 
of a universal (Nf, Nf) symmetry breaking pattern. 

When these arguments are combined with the universality of the chiral 
phase transition as advanced by Pisarski and Wilczek |5l]] in the context of 
QCD, they show that these oscillations survive the effects of matter at the 
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critical point, thereby generalizing the arguments in the vacuum. Although 
our arguments were restricted to effective QCD, they apply equally well to 
QCD in odd dimensions, as well as effective fermionic theories with spon- 
taneous symmetry breaking such as the Nambu-Jona-Lasinio models or the 
instanton models. For the latter spectra, we note that power counting is 
simplified on the hypertorus. These assertions can be readily checked. 

In the third part of these lectures, we have applied some of the Blue's 
function techniques to chiral random matrix models, in the non-universal 
regions where the 1/N expansion applies. The latters belong to a cohort of 
effective models of QCD in the long wavelength limit, that reproduce the 
microscopic sum rules. In the macroscopic limit they are consistent with 
power counting only at tree level, hence their limited range of consistency 
with the QCD Ward identities. They are however useful for addressing issues 
related to the thermodynamical limit, as well as mean-field regions in most 
of the phase transitions, barring the case of large Ginzburg windows [p.00[ . 
In their standard form they are schematic versions of the instanton model 
to the QCD vacuum. Unlike the latter, however, they can be extended to 
accommodate microscopic QCD in odd dimensional spaces where the issue 
of chiral symmetry is subtle, and the question of instantons obsolete. 

As effective models of QCD in the infrared, the chiral random matrix 
models can be used to illustrate the subtleties between the large volume 
limit versus the smallness or largeness of the QCD external parameters. In 
this respect, we have analyzed the U(l) problem, and the dependence of the 
vacuum partition function on the 9 angle. In both cases the large volume 
limit is subtle when light quarks are present. Needless to say, that lattice 
QCD simulations at finite 6 are not yet available, owing to the complex 
character of the fermion determinant. 

We have shown that the concept of addition law of Blue's functions finds 
a simple realization when the temperature or quark chemical potential are 
schematically inserted in the standard chiral random matrix models. These 
models are tied to QCD only in the context of mean-field universality, and 
may be given non-standard versions to accommodate for non-mean-field 
universalities. As effective models to the QCD phase transition, they can 
be regarded as crude variants of the Nambu-Jona-Lasinio model, or the 
instanton model, to which they naturally relate. We stress, however, their 
simplicity and also their user-easy access to the quark spectra. Under a 
phase change, the latters are expected to undergo structural changes near 
zero virtuality provided that the soft and hard modes decouple in the large 
volume limit. This is the case for infrared sensitive observables, including 
order parameters. In this sense, our comparison to the present lattice data 
at the transition point is striking [j7^ ]. We note, however, that the absence 
of thermal and Fermi motion in these models, makes their extension to most 
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observables suggestive at best. 

Chiral random matrix models of effective QCD are also interesting be- 
cause they tie some of the fundamental aspects of QCD to other models 
of disorder. In this context chiral random matrix models at finite chemical 
potential resemble closely models of disorder with dissipation. The analysis 
in terms of Blue's functions generalize handily to these cases as well. 

We have presented models with strong dissipation, currently in use to 
model quantum chaotic scattering with many open channels, as well as 
models of weak dissipation, such as the model recently discussed by Hatano 
and Nelson. The latter is disordered quantum mechanics in d dimensions, 
that is still amenable to random matrix techniques provided that freeness 
is used. We have shown that this approximation is identical to the CPA 
approximation commonly in use for hermitian disordered problems {e.g. 
Anderson model). 

The models of strong nonhermiticity discussed here, apply to a large 
number of physical phenomena. They show up in nuclear, atomic and molec- 
ular spectroscopy, mesoscopic ballistic devices, microwave cavities and the 
physics of condensed state chemical reactions, to cite a few [88]. Our method 
should shed more light on the already existing analyses. 

The current models of weak nonhermiticity offer new field of investi- 
gations that are worth pursuing with the present methods. The model of 
Hatano and Nelson [B2[ (nonhermitian analogue of the Anderson model JTJ) 
for the description of flux depinning in type II superconductors offers a good 
starting point. It is a canonical example of nonhermitian localization, that 
may become relevant to a variety of other models with a directed diffusion, 
such as: models for the growth of bacterial populations |101| ], models for 
spreading genetic mutations, models of chemical reactions, neural networks, 
turbulence, and models of growth of monetary capital. We are confident that 
Voiculescu's idea of freeness will open a vested interest in this rich array of 
problems. 
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